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magnetic  fields  radiated  during  lighting  discharges  to  the  currents  that  produce  those  fields.  A 
principal  prediction  of  this  model  is  that  the  distant  (radiation)  fields  are  directly  proportional  to  the 
current  propagating  along  the  lightning  channel,  multiplied  by  the  velocity  of  propagation.  This  paper 
examines  the  derivation  of  this  relationship  and  its  implications  in  some  detail.  We  show  that  the 
formulas  commonly  used  to  describe  the  transmission-line  model  cannot  be  correctly  applied  to  many 
lightning  processes.  A  correction  factor  is  required  that  is  significant  when  the  channel  is  not  oriented 
perpendicular  to  the  line  of  sight  to  the  observer,  unless  the  propagation  velocity  is  only  a  small 
fraction  of  the  speed  of  light.  An  important  implication  of  this  more  general  version  of  the  model  is  that 
the  current  and  velocity  decouple  so  that  they  could  be  separately  determined  by  suitable,  multistation 
measurements. 
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The  “transmission-line"  model  is  an  approximation  that  is  frequently  used  to  relate  the  electric  and 
magnetic  fields  radiated  during  lightning  discharges  to  the  currents  that  produce  those  fields.  A 
principal  prediction  of  this  model  is  that  the  distant  (radiation)  fields  are  directly  proportional  to  the 
current  propagating  along  the  lightning  channel,  multiplied  by  the  velocity  of  propagation.  This  paper 
examines  the  derivation  of  this  relationship  and  its  implications  in  some  detail.  We  show  that  the 
formulas  commonly  used  to  describe  the  transmission-line  model  cannot  be  correctly  applied  to  many 
lightning  processes.  A  correction  factor  is  required  that  is  significant  when  the  channel  is  not  oriented 
perpendicular  to  the  line  of  sight  to  the  observer,  unless  the  propagation  velocity  is  only  a  small 
fraction  of  the  speed  of  light.  An  important  implication  of  this  more  general  version  of  the  model  is  that 
the  current  and  velocity  decouple  so  that  they  could  be  separately  determined  by  suitable,  multistation 
measurements. 


I.  Introduction 

In  the  “transmission-line”  model  for  lightning  discharges 
it  is  assumed  that  the  current  propagates  along  the  channel  at 
constant  velocity  and  without  changing  shape.  That  is,  it  is 
assumed  that  the  current  advances  up  the  channel  as  a 
traveling  wave,  i(r' ,  t)  =  /(/  -  1  •  ttv),  where  v  is  the 
velocity  of  propagation.  1  is  a  unit  vector  in  the  direction  of 
the  propagation  along  the  channel  (assumed  to  be  straight), 
and  f  is  the  position  vector  of  a  point  on  the  channel.  That 
current  should  propagate  along  the  channel  was  originally 
suggested  by  Dennis  and  Pierce  [1964]  as  a  means  of 
addressing  physical  inconsistencies  in  the  Bruce-Golde 
model  [ Bruce  and  Golde,  1941],  and  this  concept  was 
developed  in  detail  by  Uinan  and  colleagues  [Uman  and 
McLain,  1969,  1970;  McLain  and  Uman ,  1971;  Uman  et  at., 
1975a,  h\  Lin  et  al.,  1980;  and  others]. 

The  term  transmission-line  model  is  attributable  to  Uman 
and  McLain  [1969,  1970],  who  also  were  the  first  to  show 
that  when  the  observation  point  is  far  from  the  channel  the 
transmission-line  model  indicates  that  the  radiated  magnetic 
field  is  proportional  to  the  current  waveform,  lit),  propagat¬ 
ing  up  the  channel  [Uman  and  McLain,  1970],  The  formula 
obtained  by  Uman  and  McLain  [1970,  equations  (11)  and 
*  (12)]  provides  a  means  of  determining  the  current  remotely 

from  measurements  of  the  radiated  magnetic  field  and  an 
*•  assumed  value  of  the  propagation  velocity.  The  equivalent 
formula  for  the  electric  field  was  applied  recently  to  return 
strokes  in  triggered  lightning,  in  experiments  in  which  cur¬ 
rent,  electric  field,  and  velocity  all  were  measured,  to  test 
the  model  [Willett  et  al.,  1988,  1989], 

An  appealing  attribute  of  the  transmission-line  model  is 
that  it  yields  a  simple,  analytic  formula  relating  channel 
current  to  radiation  field.  The  version  of  the  model  described 
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above,  however,  applies  only  to  a  straight,  vertical  channel 
and  to  an  observer  on  the  ground  fat  from  the  channel  base. 
The  current  waveform  and  propagation  velocity  must  remain 
unchanged  along  the  entire  length  of  the  channel.  Most  of 
these  restrictions  can  be  lifted  if  the  mo^el  is  applied  to 
piecewise  linear  segments  of  the  channel  [Le  Vine  and 
Meneghini.  1978a,  b].  This  approach  yields  a  computation¬ 
ally  simple  means  of  accurately  calculating  the  radiation 
fields  from  tortuous  channels  on  which  the  current  waveform 
may  change  its  magnitude,  shape,  and  propagation  velocity 
from  segment  to  segment.  (Of  course,  the  induction  and 
electrostatic  components  of  the  fields,  which  become  impor¬ 
tant  if  the  range  is  not  large  compared  to  the  overall 
dimensions  of  the  channel,  are  still  not  represented.) 

Recently,  it  has  been  shown  that  the  original  analysis  of 
the  transmission-line  model  by  Uman  and  McLain  [1970] 
was  not  completely  correct.  Rubinstein  and  Uman  [1990] 
demonstrated  that  the  “turn-on"  term  defined  by  Uman  and 
McLain  [1970,  equation  (5)]  is  incomplete.  Rubinstein  and 
Uman  [1990,  equation  (12)]  derived  a  correction  factor  that 
gives  the  exact  expression  for  the  turn-on  term.  The  original 
formula  is  an  approximation  that  applies  if  the  channel  is 
oriented  perpendicular  to  the  line  of  sight  to  the  observer  or 
if  the  propagation  velocity  is  small  compared  to  the  speed  of 
light.  This  approximation  is  the  result  of  the  assumption  that 
the  retarded  time  (/  -  R/c,  where  R  is  the  distance  from  the 
source  point  to  the  observer  and  c  is  the  speed  of  light)  is 
constant  along  the  channel  segment  [Uman  and  McLain, 
1970,  discussion  below  equation  (3)).  In  effect,  Rubinstein 
and  Uman  [1990]  relaxed  this  assumption  for  the  turn-on 
term.  However,  they  said  little  about  the  implications  of 
their  correction  to  the  application  of  the  transmission-line 
formula  [Uman  and  McLain,  1970,  equations  (1 1)  and  (12)]. 

In  particular,  the  assumption  of  constant  retarded  time 
requires  that  the  channel  segment  be  perpendicular  to  the 
line  of  sight  for  most  of  the  interesting  cases  ie.g.,  where  the 
current  waveform  has  a  steep  risetime  and  the  propagation 
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velocity  is  high).  However,  the  original  transmission-line 
formula  [Uman  and  McLain,  1970,  equations  (11)  and  (12)] 
suggests  otherwise  because  it  contains  a  geometric  factor 
representing  the  orientation  of  the  channel  relative  to  the 
observer.  Rubinstein  and  Uman  [1990]  reinforce  this  erro¬ 
neous  impression  by  discussing  only  the  turn-on  term  and 
not  the  implications  to  applications  of  the  transmission-line 
formula  itself  [e.g.,  Willett  et  al.,  1988]. 

In  the  present  paper  we  examine  the  effects  of  channel 
orientation  on  the  radiation  fields  from  a  propagating  current 
pulse.  We  show  that  the  transmission-line  formulas  can  be 
derived  in  the  time  domain  without  resorting  to  the  definition 
of  a  turn-on  term.  This  analysis  leads  to  a  correction  factor 
essentially  identical  to  that  derived  by  Rubinstein  and  Uman 
[1990]  and  also  in  agreement  with  formulas  for  the  transmis¬ 
sion-line  model  derived  by  Le  Vine  and  Meneghini  [1978a, 
b ]  based  on  a  frequency  domain  analysis.  This  correction 
factor  can  be  significant  for  channels  that  are  not  perpendic¬ 
ular  to  the  line  of  sight  to  the  observer,  such  as  might  occur 
in  models  for  the  stepped  leader,  branches  of  a  return  stroke, 
or  intracloud  portions  of  the  discharge. 

In  the  sections  to  follow,  a  derivation  will  be  given  of  the 
expressions  for  the  radiated  electric  and  magnetic  fields  from 
arbitrarily  oriented  channel  segments.  (The  analysis  in  the 
body  of  the  paper  applies  to  an  isolated  channel  segment  in 
an  unbounded  region.  The  equations  for  an  arbitrarily  ori¬ 
ented  segment  above  a  conducting  plane  are  derived  in 
Appendix  B.)  This  will  be  followed  by  a  discussion  of  the 
effects  of  the  correction  factor.  It  will  be  shown  that  this 
factor  can  resolve  some  of  the  direction  ambiguity  inherent 
in  the  original  transmission-line  formula  of  Uman  and 
McLain  [1970,  equations  (1 1)  and  (12)],  suggesting  that  very 
interesting  results  might  be  obtained  from  multistation  mea¬ 
surements  like  those  proposed  by  Thomson  [1988], 


2.  Derivation 

In  this  section  a  solution  is  presented  for  the  electromag¬ 
netic  fields  radiated  from  a  short,  arbitrarily  oriented, 
straight  channel  segment  in  the  case  when  the  current  on  the 
channel  is  the  traveling  wave,  l(t  -  1  •  r'/v).  The  geometry 
is  illustrated  in  Figure  1a.  The  formal  solution  is  a  classic 
problem  which,  in  the  context  of  the  application  to  lightning, 
has  been  developed  in  the  time  domain  by  Uman  and 
colleagues  [Uman  and  McLain ,  1970;  Uman  et  al.,  1975 b; 
Master  and  Uman ,  1983  (see  correction  by  Rubinstein  and 
Uman  [1991]);  and  others]  and  in  the  frequency  domain  by 
Le  Vine  and  Meneghini  [1978a,  b].  For  example,  the  solu¬ 
tion  may  be  written  in  the  time  domain,  using  vector  notation 
to  be  general,  as  follows  [Le  Vine  and  Meneghini,  1983]; 

Mo  f  ...  ds' 

E(r,  t)  =  [/]{/  -  (/  •  XR)XR]  — 

4?r  _ _  R 


Mo  c 
4  IT 
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segment 


[/]{/ 


3  (l-VR)XR} 


ds ' 

7F 


Mo  c2 
4tt 


[/]{/- 3 (}-XR)XR) 


ds' 

W7' 
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Fig.  la.  Geometry  used  in  deriving  the  transmission-line  formu¬ 
las.  The  vector  f  denotes  a  point  on  the  channel  segment,  r  is  the 
position  vector  locating  the  observer,  and  R(r,  r‘ )  is  the  distance 
between  these  two  points.  The  VR  is  a  unit  vector  along  Rir.  r") 
pointing  from  the  source  point  toward  the  observer. 


B(r,  t)  = 


Mo  f 

4lr  Js 


segment 


/  x  XR 
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ds' . 
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Notice  that  the  integrals  in  equations  (1)  and  (2)  are  line 
integrals  evaluated  along  the  channel  segment.  The  integra¬ 
tions  are  done  in  the  “primed”  coordinate  system,  and  in  the 
notation  employed  here,  ds'  denotes  a  differentia!  'ength 
along  the  segment  and  XR  is  a  unit  vector  pointing  from  the 
source  point  at  f  to  the  observer  at  r  (i.e.,  along  R(r,  r' )  in 
Figure  la).  The  brackets  ([  ])  in  (1)  and  (2)  denote  the 
“retarded”  value  of  the  argument  (e.g.,  see  equation  (5a)), 
and  the  dot  above  the  brackets  denotes  a  derivative  with 
respect  to  the  argument  of  the  function  enclosed  by  the 
brackets. 

Equations  (1)  and  (2)  represent  the  solution  for  the  electric 
and  magnetic  fields  due  to  a  transmission-line  current  pulse 
propagating  in  the  I  direction  (which  is  also  taken  as  the 
positive  direction  for  current  flow).  They  apply  to  an  isolated 
channel  segment  (i.e.,  the  effect  of  a  ground  plane  to 
represent  the  earth  has  not  been  included).  The  solution  is 
equivalent  to  that  obtained  by  Uman  and  colleagues.  For 
example,  (1)  and  (2)  can  be  reduced  to  the  solution  in 
spherical  coordinates  for  an  infinitesimal  dipole  presented  by 
Uman  et  al.  [1975b,  equations  (A6)  and  (A9)  in  the  appendix] 
by  recognizing  that  when  the  current  propagates  in  the  z 
direction  (i.e.,  /  =  z),  as  assumed  in  the  paper  cited  above, 
one  has  (/  •  XR)XR  =  cos  (6)fc.  I  -  (t  ■  XR)XR  =  -sin  iff) 
d,  and  /  x  XR  =  sin  (0)  d>,  where  R,  6,  4>  are  unit  vectors 
along  the  coordinate  axes  in  a  spherical  coordinate  system. 
An  alternative  derivation  of  (1)  and  (2),  showing  that  con¬ 
servation  of  charge  is  maintained  at  the  channel  end  points, 
has  been  presented  by  Le  Vine  and  Meneghini  [1983]. 

Unfortunately,  even  in  the  simple  case  of  the  transmis¬ 
sion-line  model,  general  solutions  for  the  integrals  in  (1)  and 
(2)  do  not  exist,  and  it  is  necessary  to  make  approximations. 
A  common  approximation,  which  is  useful  in  the  context  of 
radiation  from  lightning,  is  to  assume  that  the  observer  is  far 
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from  the  channel  compared  to  its  length  (i.e.,  Rif,  r‘)  »  L). 
This  is  one  of  the  approximations  made  by  Uman  and 
McLain  [1970)  in  their  original  analysis  of  the  stepped 
leader,  and  it  can  be  extended  to  longer  channels  and  to 
tortuous  channels  by  adopting  a  piecewise  linear  model  for 
the  channel  geometry  [e.g.,  Le  Vine  and  Meneghini,  1978a]. 
IfR  is  sufficiently  large,  then  only  the  terms  of  lowest  order 
in  l/R  need  be  included  in  (1)  and  (2),  in  which  case  one  has 

Mo  f  .  ,  ds‘ 

E(f,r)=-—  [/]{/- (i-VR)VR}  —  (3) 

4rr  J  segment  R 


B(r,  t)  = 


Mo  C 
4nc  J  segment 
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/  x  VR 


R 


■is' 


(4) 


where 


[/]  =  1(1  -  1  •  r'lv  -  R(r ,  r')/c).  (5a) 

Now  noting  that  ( dlds' )  =  i /  •  V’>  and  that  VR  =  -V'R, 
where  V'  is  the  gradient  operator  operating  on  the  primed 
coordinates,  one  obtains 


•  {/(/  -  /„)  -  Hi  -  th)}  (la) 


B(r.  t) 


Mo!’  Z  x  VR0 

4ttR0c  [1  -  (vlc)(l  ■  VR0)] 


•  {I(t-  tj-  I(t-  th)\  (lb) 

where  ta  =  Z  •  r'Jv  +  R(r,  r'a)!c  and  th  =  Z  •  r'h/v  +  Rir, 
r'b)lc  and  r'a  h  represent  the  ends  of  the  channel  segment. 
Notice  that  the  first  term  in  ta  h  represents  the  time  at  which 
the  current  pulse  reaches  the  indicated  end  point  of  the 
segment,  and  the  second  term  represents  the  time  required 
for  the  radiation  to  propagate  from  that  end  point  to  the 
observer.  Changes  in  R  between  end  points  have  not  been 
neglected  in  ta  or  th  because  they  can  be  important  if  the 
current  risetime  is  short,  as  will  be  shown  in  section  3. 

Now  assuming  that  i(t)  =  0  for  t  <  0,  notice  that  for  a 
short  period  after  the  radiation  initially  reaches  the  observer 
(i.e.,  t  >  tu  but  t  <  ih),  the  fields  in  (la)  and  (lb)  are 
proDortional  to  the  current  at  the  initial  end  point  (r'u)  of  the 
channel.  That  is,  in  the  time  interval,  tu<t<  th.  (la)  and 
(lb)  are 
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ds' 
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r,t  4-n-/?,)  [  1  —  (vlc)(l  •  VR0)]  ,J 


(8a) 
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Notice  that  the  effect  of  R(r,  f‘)  on  the  retarded  time 
(?  ■  VR/c  in  (5b))  is  the  term  missing  from  the  analysis  of 
Uman  and  McLain  [1970,  equation  (7)]  and  corrected  by 
Rubinstein  and  Uman  [1990,  where  it  occurs  in  a  mathemat¬ 
ical  identity  involving  the  integral  of  a  delta  function]. 
Rearranging  (5b),  we  have  the  relation 


[1  -  (vlc)0  ■  VR)]  ds 


iin 


(5c) 


which  can  be  used  in  (3)  and  (4)  to  rewrite  [/]  in  terms  of 
derivatives  with  respect  to  the  dummy  variable  of  intergra- 
tion,  s' . 

For  observation  points  sufficiently  far  from  the  channel 
(i.e.,  R(r,  r')  ?s>  L),  it  is  reasonable  to  neglect  changes  in  R 
and  VR  over  the  length  of  the  channel  segment.  Replacing  R 
and  VR  by  their  values  at  the  center  of  the  segment  (denoted 
below  by  the  subscript  “0”)  everywhere  except  in  [/]  (see 
below)  and  then  factoring  them  out  of  the  integral,  one 
obtains 
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Finally,  performing  the  integrations  indicated  in  (6a)  and 
(66),  one  obtains 
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This  is  the  important  result  noticed  by  Uman  and  McLain 
[1970]  that  generated  interest  in  estimating  the  current  in 
lightning  discharges  indirectly  from  observation  of  the  radi¬ 
ated  electromagnetic  fields.  (Notice,  however,  that  the  pro¬ 
portionality  constant  involves  the  velocity  of  propagation,  v, 
distance,  R0,  and  orientation  of  the  segment,  1,  all  of  which 
must  also  be  determined  in  order  to  obtain  a  valid  measure¬ 
ment  of  l(t).) 

Equations  (lb)  and  (86)  for  a  single  channel  segment  in  a 
boundary-free  region  correspond  to  the  transmission-line 
formulas  obtained  by  Uman  and  McLain  [1970,  equations 
(11)  and  (12))  for  a  vertically  oriented  channel  and  its  image 
in  a  conducting  plane.  Aside  from  the  generalization  made 
here  to  include  arbitrarily  oriented  channel  segments,  how¬ 
ever,  (76)  and  (86)  differ  from  the  results  presented  by 
Uman  and  McLain  [1970]  in  two  important  respects.  The 
first  is  the  presence  of  the  multiplicative  factor  [1  -  (vlc)0  • 
VR)]-1,  which  is  not  present  in  the  expressions  derived  by 
Uman  and  McLain  [1970],  The  second  is  the  dependence  of 
the  time  delays,  ta  h,  on  the  time  required  for  the  electro¬ 
magnetic  radiation  to  propagate  from  the  ends  of  the  seg¬ 
ment  to  the  observer  (i.e.,  the  terms  R(r,  r'„)/c  and  R(r, 
r’b)lc  in  the  expressions  for  ta  and  th).  These  terms  are  not 
present  in  the  expressions  presented  by  Uman  and  McLain 
[1970;  see  equations  (11)  and  (12)  and  the  discussion  below 
equation  (10)  where  “t,,”  is  defined).  Both  of  these  differ¬ 
ences  result  from  our  relaxation  of  the  assumption  of  con¬ 
stant  retarded  time. 

In  order  to  gain  insight  into  the  nature  of  (7)  and  (8)  and  to 
better  compare  with  the  transmission-line  formula  obtained 
by  Uman  and  McLain  [1970],  equations  (7)  and  (8)  will  be 
examined  in  section  3  in  the  special  case  of  a  vertical  channel 
segment.  It  will  be  shown  in  this  example  that  the  differences 
mentioned  above  can  be  quite  significant. 
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Fig.  16.  Geometry  used  to  compute  an  approximation  for  the 
time  delays  ta  and  tb  that  is  valid  when  Ra  »  L.  This  figure  also 
defines  the  spherical  coordinate  system  introduced  in  section  3. 


3.  Example 


Vertical  Channel  Segment 


In  order  to  more  readily  visualize  (7a)  and  (76)  above, 
consider  the  solutions  expressed  in  a  spherical  coordinate 
system  with  its  origin  at  the  center  of  the  channel  segment 
and  oriented  so  that  the  current  propagates  in  the  direction 
of  the  positive  z  axis,  as  illustrated  in  Figure  16.  The 
observer  is  located  at  the  point  R,  0,  <f>.  (Since  no  boundaries 
are  present,  this  represents  a  completely  general  case.  The 
solutions  for  an  arbitrarily  oriented  channel  segment  above  a 
conducting  plane  are  worked  out  in  Appendix  B.)  In  this 
case,  one  finds  that  VR0  =  R  and  1  =  z  =  cos  (9)  A  -  sin 
(0)0,  and  it  is  easily  shown  that  1  -  (1  •  V2?0)V/?0  =  -sin  (0)0 
and  that  l  x  VR0  =  sin  (0)<t>,  where  ft,  0,  and  4>  are  the 
orthogonal  unit  vectors  in  the  new  coordinate  system.  Using 
these  expressions  in  (7a)  and  (76),  one  obtains 


\Vi/£o 

E(C  t )  =  (vie) 

4-rrRn 


sin  ( 0)6 


[1  -  (vie)  cos  (0)] 


B(r ,  t)  = 


Mo 
4  ttR, 


(vie) 


■  [I(t  -  ta)  -  l(t  -  f*)]  (9a) 

sin  (0)<J> 


[1  -  (vie)  cos  (0)] 


■  [/(/  -  ta)  -  I(t  -  t„)]  (96) 


Notice  that  the  fields  in  (9a)  and  (96)  are  orthogonal  to 
each  other  and  to  the  line  of  sight.  A,  from  the  channel 
segment  to  the  observer.  Also  notice  that  the  ratio  of  the 
magnitude  of  the  electric  and  magnetic  fields  is  the  speed  of 
light  in  vaccum  (i.e.,  [E  ■  0/B  •  <£>]  =  c).  These  are  properties 
that  one  expects  to  find  in  the  radiation  components  of 
electric  and  magnetic  fields. 

Both  fields  in  (9a)  and  (96)  are  proportional  to  the  ratio 


Fig.  2.  Polar  plot  of  the  absolute  value  of  the  geometric  factor, 
sin  (0)/[l  -  (v/c)  cos  (0)],  for  several  different  values  of  v/c.  The 
angle  fl  is  zero  along  the  vertical  axis  (pointing  up)  and  increases  in 
the  clockwise  direction  (see  also  Figure  16). 


sin  (0)/[l  -  (v/c)  cos  (0)]).  In  the  equivalent  result  obtained 
by  Uman  and  McLain  [1970,  equations  (11)  and  (12)],  the 
magnetic  field,  B(r,  t),  is  proportional  only  to  sin  (0),  and  the 
factor  [1  -  (v/c)  cos  (0)]~‘  is  not  present.  In  order  to 
illustrate  the  significance  of  this  difference  the  absolute  value 
of  the  geometry  factor  sin  (0)/[l  -  (v/c)  cos  (0)]  has  been 
plotted  in  Figure  2  as  a  function  of  0  for  several  values  of 
(v/c). 

Notice  that  when  v/c  «  1,  the  expression  [1  -  (v/c)  cos 
(0)]  as  1,  in  which  case  (96)  reduces  to  the  result  obtained  by 
Uman  and  McLain  [1970]  (except  for  a  factor  of  2  introduced 
because  they  included  the  effect  of  a  perfectly  conducting 
ground  plane  in  their  analysis  (see  Appendix  B)).  In  this 
case,  the  radiation  pattern  is  identical  to  the  classical  pattern 
of  radiation  from  a  dipole  [e.g.,  Jackson,  1966].  However,  as 
the  ratio  v/c  increases,  the  radiation  pattern  changes,  be¬ 
coming  peaked  and  tilted  toward  the  direction  of  propagation 
of  the  current  pulse.  This  is  an  effect  seen  in  traveling  wave 
antennas  [e.g.,  Jordan  and  Balmain,  1968]. 

Of  course,  the  magnitude  of  the  fields  in  (9a)  and  (96)  also 
depends  directly  on  the  ratio  (v/c),  as  shown  by  Uman  and 
McLain  [1970],  That  is,  not  only  does  the  radiation  tend  to 
peak  in  the  direction  of  propagation  along  the  channel 
segment  as  v/c  — *  1  but  also  the  magnitude  of  the  radiation 
itself  increases.  This  is  illustrated  in  Figure  3,  where  the 
absolute  value  of  the  quantity  (v/c)  sin  (0)/[l  -  (v/c)  cos  (0)] 
has  been  plotted  for  the  same  values  of  v/c  as  used  in 
generating  Figure  2. 

Finally,  note  that  the  radiation  patterns  in  Figures  2  and  3 
are  rotationally  symmetric  about  the  z  axis  because  in  the 
geometry  being  treated  here,  there  is  no  dependence  on  the 
coordinate,  <6. 

In  addition  to  the  explicit  dependence  of  the  fields  on  the 
segment-observer  geometry,  as  described  above,  there  is 
also  an  implicit  dependence  included  in  (9a)  and  (96)  in  the 
factor,  I(t  -  ta)  -  I(t  -  tb).  This  is  demonstrated  in 
Appendix  A,  where  it  is  shown  that  when  LIR0  <K  1,  one 
obtains 


L 

ta  =  R0/c  -  —  [1  -  (v/c)  cos  (0)] 
2  v 


(10  a) 
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Fig.  3.  Polar  plot  of  the  absolute  value  of  the  factor,  [(u/c)  sin 
(0)]/[l  -  (vie)  cos  (0)],  for  several  different  values  of  vie.  The  angle 
6  is  zero  along  the  vertical  and  increases  clockwise. 


L 

tb  =  Rq/c  +  —  [1  —  (v/c)  cos  (0)]  (106) 


and  we  define 

T  =  ‘b~  ‘a=  U-M[l  ~  (v/c)  COS  (0)]  (10c) 

In  (10a),  (106),  and  (10c),  L  is  the  length  of  the  channel 
segment. 

As  mentioned  earlier,  ia  and  th  consist  of  the  time 
required  for  the  current  pulse  to  propagate  to  the  appropriate 
end  of  the  segment,  plus  the  time  required  for  electromag¬ 
netic  radiation  to  cover  the  distance  from  the  respective  end 
point  to  the  observer.  Since,  in  general,  neither  of  these 
times  is  the  same  for  the  two  end  points,  one  obtains  the 
nonzero  net  time  difference  r  defined  in  ( 10c).  The  value  of 
t  is  the  difference  in  the  retarded  times  associated  with  the 
arrival  of  the  current  pulse  at  the  two  end  points  of  the 
channel  segment. 

Notice  that  the  first  term  in  r  (i.e. ,  Llv)  is  just  the  time 
required  for  the  current  to  propagate  from  one  end  of  the 
channel  to  the  other  end.  In  the  original  work  on  the 
transmission-line  model  only  this  term  was  taken  into  ac¬ 
count.  As  a  result,  Uman  and  McLain  [  1970,  see  equations 
(11)  and  (12)  and  the  discussion  below  equation  (10)  where  rv 
is  defined)  ignored  the  dependence  on  the  angle  0  contained 
in  the  second  term  in  (10c)  (i.e.,  -( Lie )  cos  0). 

The  effect  that  this  implicit  dependence  on  the  channel- 
observer  geometry  can  have  on  the  radiation  depends  on  the 
shape  of  the  current  waveform,  /(f).  In  particular,  when  1(1) 
has  a  rapid  risetime,  as  is  typical  of  the  current  in  return 
strokes,  small  variations  in  r  can  cause  large  changes  in  the 
observed  fields.  In  this  case,  (9)  and  (10)  can  be  significantly 
different  from  those  in  which  rv  =  Llv). 

To  illustrate  this,  the  difference  l(t  -ta)  -  l(t  -  th)  has 
been  plotted  in  Figure  4 6  for  several  different  values  of  the 
angle  0.  To  generate  these  examples,  it  has  been  assumed 
that  L  =  100  m  and  v  =  10s  m/s.  An  exponential  model  has 
been  chosen  for  the  current  waveform,  as  illustrated  in 
Figure  4a:  /(f)  =  /0  [exp  (-at)  -  exp  (~bt) J,  where  a  = 
2  x  104  s"1  and  b  =  3.5  x  106  s_l.  This  is  the  current 
waveform  suggested  by  Bruce  and  Golde  [1941]  with  param¬ 
eters  suggested  by  Dennis  and  Pierce  [1964]  to  describe 


TIME  (microseconds) 

Fig.  4.  (a)  The  current  waveform  used  in  ( b )  to  compute  the 
factor  [/(/)  -  l(t  -  r)]//0  for  various  values  of  0.  The  0-dependent 
element  of  this  expression  is  r  =  (L/v)[  1  -  (vie)  cos  (0))  (equation 
<  10c)).  The  solid  curve  (0  =  90°)  corresponds  to  the  result  obtained 
by  Uman  and  McLain  [1970]  for  all  values  of  0. 


return  strokes,  except  that  b  has  been  modified  to  reflect 
modem  measurements  that  suggest  a  more  rapid  risetime 
[e.g.,  Weidman  and  Krider,  1980,  1982]. 

Equation  10c  shows  that  r  is  a  minimum  at  0  =  0°  and 
increases  to  a  maximum  at  0  =  180°.  The  consequence  is 
evident  in  the  four  curves  in  Figure  46.  For  example,  the 
curve  for  0  =  150°  has  the  largest  peak  value  because  the 
longer  retarded  time  difference  permits  /(/  -  ta)  to  grow 
more  before  the  term  -/(/  -  tb)  arrives  and  begins  to  cancel 
it.  Similar  results  would  be  obtained  for  other  choices  of  L. 
v  and  /(f),  although  changing  the  current  waveform  would 
change  the  specific  shape  of  the  curves. 

As  described  above,  the  expression  [1  -  (vie)  cos  (0)] 
affects  the  radiated  fields  in  two  ways:  explicitly  through  the 
geometry-dependent  factor  sin  (0)/[l  -  (v/c)  cos  (0)],  as 
illustrated  in  Figures  2  and  3,  and  implicitly  through  the 
difference  I(t  -  ta)  -  I(t  -  tb),  as  illustrated  in  Figure  4. 
The  combined  effect  of  these  two  terms  is  shown  in  Figure  5, 
where  the  complete  solution  has  been  plotted  for  several 
values  of  0,  using  the  parameters  described  above.  For  this 
purpose  a  normalized  field  magnitude,  F(0,  f),  has  been 
defined  as  follows: 


F(0,  f)  =  |£(f,  f)| 


VM(/Eo 

4rrR0 


(vlc)I o 
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malized  field,  defined  in  (1 1),  has  been  plotted  after  dividing 
by  sin  (0),  a  factor  common  to  both  solutions. 

The  Uman  and  McLain  result  corresponds  to  the  solid 
curve  in  Figure  6,  and  it  applies  lor  all  0  since  sin  (9)  is  the 
only  9  dependence  in  their  solution.  Dividing  (9a)  and  (9b) 
by  sin  (9)  leaves  the  implicit  and  explicit  dependences  on  [1 
-  (vie)  cos  (0)]  added  by  the  present  analysis.  Our  result  has 
been  plotted  for  several  values  of  9 ,  using  the  parameters 
described  above.  It  is  clear  from  Figure  6  that  (9a)  and  (9b) 
can  differ  significantly  from  the  Uman  and  McLain  solution 
for  values  of  0  other  than  90°. 


Fig.  5.  Complete  solution  (equations  (9)  and  (10)),  using  the 
current  waveform  of  Figure  4a,  for  several  different  angles  6.  The 
normalized  field  F(8,  /),  defined  by  equation  (II),  has  been  plotted. 


=  m  oh 


(v/c)I0\. 


Comparing  Figures  2  and  4b ,  it  is  clear  that  the  explicit 
and  implicit  contributions  affect  the  received  fields  in  dif¬ 
ferent  ways.  The  implicit  dependence  results  in  an  increase 
in  the  amplitude  as  9  increases,  because  the  time  delay  is 
increased  (Figure  4b).  On  the  other  hand,  the  effect  of  the 
explicit  dependence  (Figure  2)  is  to  increase  the  amplitude  to 
a  peak  near  9  =  60°  and  then  to  decrease  it  with  further 
increases  in  angle.  It  is  clear  from  Figure  5  that  the  explicit 
dependence  dominates  in  this  case,  resulting  in  a  net  de¬ 
crease  in  amplitude  for  angles  greater  than  about  60°. 

Equations  (9)  and  (10)  are  equivalent  to  the  results  ob¬ 
tained  by  Uman  and  McLain  [1970,  equations  (11)  and  (12)] 
when  9  =  90°,  because  [1  -  (v/c)  cos  (0))  =  1  in  this  case. 
However,  for  other  values  of  0the  two  can  be  quite  different. 
This  is  illustrated  in  Figure  6,  where  the  two  solutions  have 
been  plotted  together.  (Uman  and  McLain’s  solution  has 
been  divided  by  2  to  compensate  for  the  effect  of  the 
conducting  plane  in  their  analysis.)  In  each  case,  the  nor- 


(mcroseconas) 

Fig.  6.  Comparison  of  the  Uman  and  McLain  [19701  solution, 
represented  by  the  solid  curve  for  all  9  and  equations  (9)  and  ( 10)  for 
several  different  angles,  using  the  current  waveform  of  Figure  4a. 
F(  6,  t)/sin  6  has  been  plotted  to  remove  the  common  factor  from 
both  normalized  fields.  The  two  solutions  are  the  same  when  8  = 
90°. 


4.  Apparent  Singularity 

In  addition  to  the  characteristics  of  (7)  and  (9)  pointed  out 
above,  notice  that  when  vie  =  1,  it  is  possible  for  the 
expression  [1  -  (v/c)(l  •  V/f0)]  in  (la)  and  (lb)  (or 
equivalently,  the  expression  [1  -  (v/c)  cos  (9)]  in  (9a)  and 
(9b))  to  be  zero.  This  happens  when  the  observer  is  along 
the  axis  of  the  channel  segment  in  the  direction  of  propaga¬ 
tion  of  the  current  waveform  (i.e.,  1  =  VR0  or  0  =  0).  At  first 
glance  it  appears  that  the  fields  are  unbounded  in  this  case, 
since  sin  (0)/[l  -  cos  (0)]  blows  up  as  0  goes  to  zero. 
However,  the  singularity  is  only  apparent,  and  both  fields 
actually  go  to  zero  along  the  axis  of  the  segment. 

To  see  this,  again  consider  a  spherical  coordinate  system 
with  its  origin  at  the  center  of  the  channel  segment  and 
oriented  such  that  !  =  z.  It  is  desired  to  examine  the  fields 
(equations  (9a)  and  (9b))  near  the  axis  (0  =  0)  in  the  special 
case  v/c  =  I.  Since  v  =  c,  one  can  always  get  close  enough 
to  the  axis  so  that  the  retarded  time  difference,  r  (equation 
(10c)),  is  small,  no  matter  how  long  the  channel  segment 
becomes  (i.e.,  whatever  the  magnitude  of  Lie).  Assuming 
that  ris  small,  and  using  ta  =  R0lc  -  r/2  and  tb  -  R0lc  + 
r/2  from  (10a),  (10 b),  and  (10c),  we  can  expand  l(t  -  ta) 
and  l(t  -  tb)  in  power  series  about  t  =  t  -  Role.  To  first 
order  in  t  this  yields 


/(/  -  tj  -  I(t  -  th)  s  T  -  /(f)  (12a) 

at 

t  =  ( 

Hence  close  to  the  axis  the  normalized  field  (equation  (11) 
becomes 

sin  (0) 

0  =  -  -  ■■■  - -  [/(/  -  tj  -  l(t  -  /„)]//„ 

LI  -  (vie)  cos  0] 


=  (L/vI0)  sin  (0)  —  /(/)  (12 b) 

at 

t  =  t 

which  goes  to  zero  as  0  — »  0.  Thus  both  E(r,  t)  and  B(r,  t) 
go  to  zero  along  the  axis  of  the  channel  segment,  and  there 
is,  in  fact,  no  singularity. 

One  can  reach  the  same  conclusion  using  the  following 
argument.  Notice  that  when  0=0  and  v/c  —  1,  the  current 
pulse  leaving  the  initial  end  of  the  channel  segment  (i.e.,  at  f 
=  r'a)  and  the  electromagnetic  field  radiated  from  this  end 
(i.e.,  the  term  corresponding  to  I(t  -  ta)  in  (9a)  and  (9b) 
propagate  toward  the  observer  along  the  same  ray  path  and 
at  the  same  velocity.  Consequently,  in  this  special  case  they 
both  reach  the  other  end  of  the  segment  at  the  same  time. 
When  the  current  pulse  reaches  this  end  of  the  segment,  a 
second  pulse  of  radiation  is  emitted  (the  term  associated  with 
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Hr  -  tb)  in  (9 a)  and  (96).  The  radiation  from  this  end  is 
equal  in  amplitude  but  opposite  in  sign  to  the  radiation 
arriving  from  the  initial  end  and  cancels  the  arriving  radia¬ 
tion.  Hence  the  net  electromagnetic  field  radiated  to  the 
observer  is  zero.  In  effect,  the  two  singularities  in  (9a)  and 
(96)  cancel  each  other,  leaving  the  fields  bounded. 


5.  Very  Short  Segment 

The  radiated  fields  (equations  (9a)  and  (96))  reduce  to  an 
interesting  special  form  in  the  case  of  a  very  short  channel 
segment,  as  might  be  appropriate,  for  example,  in  modeling 
radiation  from  an  individual  step  in  the  stepped  leader 
process.  Assume  that  the  segment  is  sufficiently  short  that 
L/v  «  Tr ,  where  r,  is  the  risetime  of  the  current  pulse. 
From  (10c)  the  retarded  time  difference  T«  Tr,  because  [1 
-  ( v/c )  cos  0]  <  2  for  any  v  <  c.  Now  expand  /(/)  in  a  power 
series  in  r,  as  done  in  ( (2a)  above,  keeping  only  terms  of 
lowest  order  in  r.  Inserting  (126)  into  (9 a)  and  (96),  one 
obtains 


E(r,  t)  =  0 


f  — . 

vmV£o 

4itR0 


d 

( Lie )  sin  (6)  —  /(/) 
dt 


(13a) 


B(r,  t)  =  4>  (He)  sin  (6)  —  l(t) 

4-ttRq  dt 


(136) 


As  one  would  expect  in  this  limit,  the  radiation  pattern  is 
identical  to  that  of  a  classical  dipole  [ Jackson ,  1966;  Uman  et 
al.,  19756,  equations  (A6)  and  (A8)],  regardless  of  the  ratio 
vie.  Also  notice  that  the  radiation  is  directly  proportional  to 
the  derivative  of  the  current  waveform,  rather  than  being 
proportional  to  the  current  waveform  itself  as  is  the  case  (at 
least  initially)  with  a  longer  segment  (i.e,  (7)  or  (9)).  Finally, 
notice  that  in  this  case,  the  magnitude  of  the  fields  is  no 
longer  proportional  to  the  velocity  of  propagation  of  the 
current,  v,  but  is  now  diiectly  proportional  to  the  length,  L, 
of  the  channel  segment. 

Also  notice  that  for  this  very  short  segment  the  depen¬ 
dence  on  [  1  -  (v/c)  cos  (8)]  has  entirely  disappeared.  As  the 
channel  element  gets  shorter  and  shorter,  the  explicit  and 
implicit  dependences  on  this  expression  become  nearly 
equal  and  tend  to  cancel  each  other,  indicating  both  both 
dependences  are  needed  for  a  self-consistent  solution.  This 
cancellation  is  illustrated  in  Figure  7,  where  the  normalized 
fields  have  been  plotted  (using  (9)  and  (1 D)  for  several  values 
of  8  and  for  decreasing  values  of  the  element  length,  L.  (As 
in  Figure  6,  F(8,  t)l  sin  (8)  has  been  plotted  to  emphasize  the 
effects  of  the  expression  [1  -  (v/c)  cos  (0)].)  Notice  the  clear 
dependence  on  0  for  L  =  20  m  and  the  obvious  lack  of 
dependence  on  this  factor  when  L  =  0.8  m.  It  is  also 
apparent  from  Figure  7  that  the  peak  amplitude  attained  by 
the  fields  is  decreasing  in  proportion  to  L. 

Finally,  notice  that  in  order  for  (13a)  and  (136)  to  be 
applicable  to  lightning,  it  is  not  sufficient  that  the  channel  be 
short.  In  particular,  in  order  to  neglect  higher  order  terms  in 
the  power  series  expansion,  it  was  necessary  to  assume  that 
L/v  was  small  compared  to  rr.  Thus  at  least  in  the  context  of 
a  transmission-line  model,  use  of  the  dipole  approximation 
requires  that  the  time  required  for  the  current  pulse  to 
propagate  the  length  of  the  channel  segment  be  small  com¬ 
pared  to  the  risetime  of  the  current  waveform  itself.  This  is 


Fig.  7.  Effect  of  channel  length,  L,  on  the  radiation  from  a 
straight  segment.  For  each  value  ofZ.,  F(6,  t)/ sin  0has  been  plotted 
(based  on  equations  (9),  (10),  and  (1 1))  for  several  different  angles,  8, 
using  the  current  waveform  of  Figure  4a.  Note  that  Figure  6  shows 
the  same  function  for  L  =  100  m. 


in  agreement  with  the  earlier  result  of  Krider  and  Radda 
[1975], 


6.  Consequences  for  Remote  Sensing 

Recently,  an  experiment  has  been  proposed  by  Thomson 
[1988]  to  use  wideband,  multistation  measurements  of  the 
radiated  fields  to  determine  the  location,  channel  orienta¬ 
tion,  current  waveform,  and  velocity-current  product  for 
certain  processes  in  the  lightning  discharge.  It  has  generally 
been  assumed  that  the  radiation  fields  only  contain  informa¬ 
tion  about  the  magnitude  of  the  velocity-current  product  and 
no  information  on  the  direction  of  propagation  of  the  current 
waveform  along  any  given  channel  segment.  This  is  because 
in  the  Uman  and  McLain  [1970,  equations  (11)  and  (12)] 
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TABLE  1.  Ambiguity  in  Measurement  of  Electric  Field 


Direction  of 
Propagation 

Current 

Sign 

of 

£.v«> 

Ii 

Factor 

-  iv/c)(1  •  ] 

Sign 

Direction 

r 

T 

- 

U 

-  (v/c)  cos  (0)] 

T 

- 

l 

+ 

ii 

-  (vie)  COS  (0)1 

i 

l 

+ 

U 

+  (vie)  cos  I0)j 

i 

~ 

r 

H 

+  (v/c)  cos  <0)1 

version  of  the  transmission-line  formulas,  current  and  prop¬ 
agation  velocity  appear  only  as  a  product.  However,  in  the 
formulas  derived  in  the  present  paper,  the  correction  factor 
[I  -  ( v/c)0  ■  V/?)} ~ 1  resolves  this  ambiguity  and  provides 
the  possibility  that  current  and  velocity  might  be  individually 
determined  with  suitable  multistation  measurements. 

The  nature  of  the  ambiguity  is  illustrated  in  Table  1  for  the 
special  case  of  a  vertical  channel  over  a  perfectly  conducting 
ground  plane  (see  Appendix  B).  The  first  column  in  the  table 
indicates  the  direction  of  propagation  of  the  current  wave¬ 
form  on  this  channel.  This  is  the  direction  of  the  unit  vector, 
/.  (The  velocity,  v,  is  always  a  positive  number  in  the 
formulation  in  this  paper.)  The  second  column  in  Table  1 
indicates  the  sign  of  the  current  flowing  in  the  direction 
indicated  in  the  first  column.  A  plus  sign  means  positive 
charge  moving  in  the  direction  of  the  arrow,  and  a  minus 
means  negative  charge  moving  in  this  direction.  The  third 
column  indicates  the  direction  of  the  current  vector,  iff*,  r) 
=  ll(r' ,  r),  determined  by  the  entries  in  the  first  two 
columns. 

The  fourth  column  in  Table  1  gives  the  sign  of  the  vertical 
component  of  electric  field  measured  on  the  surface,  given 
the  conditions  indicated  in  the  first  two  columns.  Here  a  plus 
sign  means  that  the  electric  field  points  vertically  upward, 
away  from  the  surface  (as  though  there  were  negative 
charges  aloft),  and  a  minus  means  that  the  electric  field 
points  downward,  toward  the  surface.  The  signs  indicated  in 
this  column  are  obtained  directly  from  equation  'BM. 

Notice  that  each  sign  in  the  fourth  column  of  Table  1 
appears  twice.  For  example,  positive  current  propagating 
upward  and  negative  current  propagating  downward  both 
yield  negative  electric  fields  on  the  surface.  Throughout  the 
table  the  direction  of  the  electric  field  is  opposite  to  the 
direction  of  the  current  vector,  indicated  by  the  arrows  in  the 
third  column,  as  one  would  expect  from  (1).  Thus  it  is  clear 
that  one  cannot,  based  on  the  sign  of  the  electric  field  alone, 
determine  the  direction  of  propagation  of  the  current  wave¬ 
form. 

In  the  original  version  of  the  transmission-line  formulas 
[e.g.,  Uman  and  McLain,  1970;  Thomson,  1988],  cases  with 
electric  fields  of  the  same  sign  in  our  table  also  have  the  same 
magnitude,  and  therefore  are  indistinguishable.  However,  in 
the  solution  presented  here,  the  correction  factor  [1  - 
( v/c){l  ■  VR0)]~*  introduces  an  additional  dependence  on 
the  direction  of  propagation,  /,  that  resolves  this  ambiguity. 
The  form  that  this  correction  factor  takes  in  each  of  the  four 
cases  is  indicated  in  the  last  column  of  Table  1 . 

Notice  that  the  correction  factor  is  different  in  the  previ¬ 
ously  ambiguous  cases,  because  the  sign  of  /  •  rRa  is 
different.  Thus  the  amplitude  of  the  electric  field  will,  in 
general,  be  different  in  these  cases.  Since  this  amplitude 
difference  depends  on  the  angle  0,  it  may  be  possible,  by 


making  measurements  at  several  locations  having  different 
values  of  cos  (0),  to  distinguish  upward  and  downward 
propagation  of  the  current  waveform.  Furthermore,  since 
the  correction  factor  also  depends  on  the  velocity,  r,  it  may 
even  be  possible  to  separately  determine  the  velocity  of 
propagation  and  the  magnitude  of  the  current.  Such  a 
possibility  would  be  of  great  interest,  for  example,  in  the 
case  of  leader  steps,  for  which  neither  the  velocity  nor  the 
direction  of  propagation  is  known  [e.g.,  Krider  and  Radda, 
1975], 

7.  Discussion 

Equations  (7 a)  and  (lb)  are  approximate  expressions  that 
apply  when  the  observation  point  is  far  from  the  channel 
segment.  However,  the  presence  of  the  factor  [1  -  (vlc)U  ■ 
Vf?0)]_l  does  not  appear  to  be  a  consequence  of  this 
approximation.  In  particular,  it  is  possible  to  obtain  an  exact 
solution  (i.e . .  no  approximations  needed  to  evaluate  the 
integrals  in  (I)  and  (2))  in  the  special  case  v  =  c  [Le  Vine  and 
Meneghini,  19786],  and  the  formulas  obtained  in  this  case, 
contain  the  factor  [1  -  (/  •  V/?0)]~'.  This  latter  solution  is 
valid  independent  of  distance  from  the  channel;  conse¬ 
quently,  the  factor  [1  -  (/  •  V /? 0 ) ]  “ 1  is  present  in  the 
radiation  component  of  the  fields  at  all  distances  R(r,  r'„) 
even  very  close  to  the  channel.  Furthermore,  the  exact 
solution  does  not  appear  to  be  singular  in  v/c  because  letting 
R(r/r'Q )  — *  *  in  this  solution  yields  (la)  and  (lb)  with  v 
replaced  by  c  [Le  Vine  and  Meneghini,  1 978/?]. 

Finally,  notice  that  the  factor  [1  -  (vlc)\l  -VR0))  1  is  not 
unusual,  appearing  in  other  applications  of  Maxwell's  equa¬ 
tions.  For  example,  this  factor  is  present  in  the  theory  for 
traveling  wave  antennas  [e.g.,  Jordan  and  Balmain,  1968] 
and  also  appears  in  the  expression  for  the  Lienard-Wiechert 
potentials  used  to  describe  radiation  from  a  moving  charge 
[Jones,  1964;  Meneghini,  1984J.  This  factor  also  appears  in 
the  correction  that  Rubinstein  and  Uman  [1990]  made  to  the 
original  theory  of  Uman  and  McLain  [1970],  as  part  of  their 
analysis  of  the  turn-on  term. 

In  summary,  the  conventional  formulas  used  to  relate 
current  and  radiation  fields  in  the  transmission-line  model 
must  be  modified  to  include  the  amplitude  factor  [  1  -  ( v/c)(I 
•  V/?0)]_l.  It  is  clear  from  Figure  2  that  this  factor  can  be 
significant  in  applications  to  radiation  from  lightning.  When 
v  «  c ,  this  factor  is  nearly  unity,  and  the  radiation  pattern 
reduces  to  the  classical  pattern  of  radiation  from  a  dipole  In 
this  case,  the  original  formulas  due  to  Uman  and  McLain 
[1970]  are  applicable.  This  correction  factor  is  also  likely  to 
be  very  nearly  unity  when  considering  radiation  from  near 
the  base  of  return  stroke  channels  to  an  observer  on  the 
ground,  a  case  in  which  the  channel  is  close  to  vertical, 
hence  perpendicular  to  the  line  of  sight  (i.e.,  /  •  VR(t  s  0). 
However,  for  radiation  from  elevated  points  on  the  channel 
and/or  from  channel  segments  of  arbitrary  orientation,  as 
might  be  encountered  in  modeling  radiation  from  leader 
steps,  channel  branches  or  intracloud  processes,  /  •  V/?0  f 
0  in  general,  and  the  correction  factor  [  1  -  (v/c)(l ■  V/?0 ) ]  " 1 
can  be  significant. 

Appendix  A 

Consider  a  channel  segment  of  length  L  oriented  along  the 
z  axis  of  a  spherical  coordinate  system  as  indicated  in  Figure 
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lb.  The  origin  is  at  the  center  of  the  segment.  The  purpose 
of  this  appendix  is  to  compute  the  time  delays  tu  =  Ralc  - 
I/(2r)  and  tb  =  RhU-  +  L/(2v)  in  the  special  case  R„  »  L. 
Thus  letting  0  be  the  angle  R0  makes  with  the  z  a.  .s,  one 
obtains  the  following  expressions  from  standard  trigonomet¬ 
ric  rules: 

Ru  =  /?o\ 1  +  (L/R0)  cosT»T+  (E/2/?0)2  (AD 


Rh  =  R 0\l  -  (L/R0)  cos  (0)  +  (L/2R0) '  <A2) 

Expanding  the  square  root  in  a  power  series  keeping  only  the 
first  order  terms  in  the  small  parameter  (LIR0).  one  obtains 

Ru  s  R()  +  (L/2)  cos  (it)  (A3) 

Rh  s  R{)  -  (L/2)  cos  (0)  ( A4) 

The  error  is  of  the  order  of  L://?0.  Using  (A3)  and  (A4),  it 
follows  that 

L 

ra=  ~L/(2v)  +  Ra/c  =  R0h  [1  -(i</r)  cos  (0)]  (A5) 


L 

th  =  +L(2r)  +  Rhlc  =  R0/c  +  —  [1  -  (vie)  cos  (  0  )].  ( A6) 

2r 


Appendix  B 

Since  most  applications  to  lightning  are  for  measurements 
made  on  the  surface  of  the  Earth,  equations  are  derived  in 
this  appendix  for  the  radiation  fields  seen  by  an  observer  on 
the  surface  when  the  source  is  'over  a  perfectly  conducting 
ground  plane.  This  problem  is  solved  by  applying  the  method 
of  images.  In  particular,  the  solutions  are  obtained  by  adding 
to  (I)  and  to  (2)  the  contribution  due  to  an  image  current, 
f/(rj,  t )  =  —  /( r  -  i,  ■  r)/v ),  where  r =  r'  -  2 z'z  and  /,  = 
[/  -  2 (/  •  2)21-  The  solutions  for  the  fields  due  to  the  image 
current  follow  the  procedure  outlined  in  (3)  through  (7)  in  the 
text  and  are  added  to  the  solutions  already  presented  there. 

When  the  observer  is  on  the  surface,  the  solutions  simplify 
greatly  because  on  the  surface  the  electric  field  has  only  a 
component  perpendicular  to  the  surface  and  the  magnetic 
field  has  only  a  component  tangential  to  the  surface.  The 
solutions  in  this  case  can  be  written  quite  simply  in  terms  of 
(7a)  and  (7b).  Letting  ft  denote  a  unit  vector  perpendicular 
to  the  surface  and  letting  Es(r,  t)  and  Bs(r,  r)  denote  the 
fields  on  the  surface,  one  obtains 

E,(r,  t)  =  2[h  •  E(r.  f)]n 


where  E(r,  t)  and  B(r.  t )  are  the  solutions  given  in  (7u)  and 
(lb),  respectively,  and  where  (Bl)  is  obtained  from  (7a) 
with  the  vector  identity  d  x  (h  x  c)  -  bid  ■  C)  -  ctd  ■  b). 

In  order  to  compare  with  other  forms  of  the  solution  that 
have  appeared  in  the  literature,  consider  a  spherical  coordi¬ 
nate  system  with  its  origin  at  the  center  of  the  channel 
segment,  as  in  section  3,  but  oriented  so  that  its  axis  is 
perpendicular  to  the  surface  In  -  r).  Then,  noting  that 
x  In  x  V/?01  =  -sin  (0)0,  (Bl)  and  (B2)  can  be  written 

pui'  sin  (  0  K  0  ■  I  )z 
2t tRu  [I  -  (i'/f  M /  -  V/?0i) 


•  [tu  -  tu)  -  lit  -  th Ij  (B3t 

_  _  Mo1'  [c  x  (/  x  VRU)]  X  : 

’  2-t tR0c  [ I  -  (vlc)(l  •  T/?„|] 

•  [/(/  -  t a)  -  /(/  -  /(,)].  ( B4) 

Note  that  the  angle  0  is  larger  than  90°  when  the  real  source 
is  above  the  ground  plane. 

In  addition,  if  the  channel  segment  is  oriented  parallel  to 
the  z  axis  (i.e..  /  =  ±z),  then  our  spherical  coordinate 
system  becomts  identical  to  that  of  Figure  Ife,  and  the 
results  simplify  further  (e.g..  /  •  T/?,,  =  i  cos  ((/).  [  z  x  (/  x 
V/?0|]  x  ;  =  ±  sin  (0)d>,  and  1  ■  0  =  7  sin  (0)).  One  obtains 

Mnv  sin"  (0  )z 

EAr,  /)=-<-  - — —  - - - 

2nR0  [  I  *  (i'/c)( cos  0  )] 


•  [/( /  -  ta)  -  Ht  -  th)]  *  B5 ) 
M  i)i'  sin  ( 0  )<f> 

BJr.  t)=  ±  -  r—z - 7 - r 

2t tR0c  [1  -i-  (v/c Kcos  0)] 

•  UU  -  ta)  -  lit  -  th)]  (B6) 

Here  the  upper  sign  applies  to  an  upward  propagating 
current  pulse  and  the  lower  to  downward  propagation. 
(Recall,  however,  that  the  positive  direction  for  current  flow 
is  taken  as  the  propagation  direction,  so  that  the  fields  do  not 
reverse  sign  unless  the  current  flow  reverses  direction.) 
Except  for  the  explicit  correction  factor  [  1  -  (u/cHcos  fl)|  1 
and  the  additional  ^dependence  implicit  in  tu  and  th.  (B6)  is 
identical  to  that  derived  by  Uman  and  McLain  [1970], 
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Bs (r,  t )  =  2 [fi  x  B(r,  /)]  x  h 
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